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Abstract 

We study the nucleon-nucleon interaction in the isoscalar-scalar channel 
using the chiral unitary approach. The i-matrix of the pion-pion scattering in 
this channel is summed up to all orders using the B-S equation. We find that 
the calculated results at long distances are close to those of the cr-exchange 
interaction. In addition, there appears a shorter range repulsion in this chan- 
nel. 

1 Introduction 



The intermediate range attraction in the NN interaction has been traditionally 
described by the cr-exchange in the meson exchange picture. |IJ It has also been 
noticed |2j, |], |J that box diagrams with two-pion exchange and intermediate A 
excitation lead to an intermediate range attraction. A weakened a exchange together 
with these box diagrams has also been used to describe the intermediate range NN 
attraction. ^ 

With the success of chiral perturbation theory (xPT) in the meson-meson and 
meson-baryon sectors, || attempts to extend these ideas to the NN sector have been 
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pursued. || [7|, [£| In a recent paper, |J following the line of Ref. [8], the peripheral 
NN partial waves are studied within a chiral scheme, using the meson- meson and 
the meson-baryon interaction Lagrangians. With this input, together with more 
conventional processes with A box diagrams, p exchange and other meson exchange, 
a good reproduction of the NN data for L > 2 partial waves is obtained. 

A striking feature of Ref. [9] is that the exchange of two interacting pions in 
the scalar- isoscalar channel (the a channel) leads to a repulsion (although weak), 
instead of the commonly accepted attraction from a exchange. Simultaneously, from 
the box diagrams with intermediate A, an attraction is obtained with the range and 
strength of the standard a exchange of the boson exchange models. Two basic 
approximations lead to these results. First, the nir isoscalar interaction is used only 
to lowest order in ^PT. Second, no form factors are considered for the A box 
diagrams. 

In the present work we wish to reconsider this idea and go further by treating 
the exchange of two interacting pions in the isoscalar channel in a nonperturbative 
way. The approach followed here for the inr interaction produces a o pole in the 
complex plane in the physical region s > 4m^. Then the analytical extrapolation of 
the model is used in order to find the strength of the isoscalar exchange interaction 
for the situation s < 0, which one encounters in the iViV scattering problem. 

The a meson has been rather problematic, with ups and downs in the particle 



data tables, where it has been once again welcome. |T0[ Some analyses of the ixix 



data rely on the a pole in the Tin t matrix. [11, 12, O Theoretical models for the 



7r7r interaction based on meson exchange fllj] find also a pole in the t matrix for the 
a meson. Yet, of relevance to the present work is the fact that the a pole is found in 
recent chiral nonpertubative approaches which have developed independently. 



T6| In Ref. [15] the inverse amplitude method using the lowest order and second order 



chiral Lagrangians is formulated and it leads to a pole position for a at 440 — i245 
MeV. In Ref. [16] the Bethe Salpeter equation, with a cutoff in the loops fitted to 
the data of the scalar sector, is shown to reproduce accurately the data in the scalar 
sector up to around = 1.2 GeV, using only the lowest order chiral Lagrangian 
as input. In this case, a appears as a pole at 469 — 2203 MeV [Ref. [16], erratum]. 
The approach of Ref. [16], using coupled channels, is able to reproduce the a and 
fo (980) resonance in the scalar-isoscalar sector and the ao(980) resonance at L = 0, 
1=1. The approach of Ref. [15] using only the irn channel could only generate a in 
the scalar channel, but p and K* in the vector channel were accurately reproduced. 

A generalization of the inverse amplitude method of Ref. [15], incorporating 
coupled channels as in Ref. [16] and using the lowest and second order chiral La- 
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grangians, is given in Ref. [17], and there, all meson-meson data up to yfs = 1.2 
GeV are accurately reproduced. In the latter work, the a pole appears at 442 — ?225 
MeV. 

A different chiral approach is employed in Ref. [18], using the N/D method. 
This explicitly intoduces resonances on top of the lowest order chiral Lagrangian. 
That approach, which includes contributions from the physical and unphysical cuts 
through dispersion relations, produces good results up to y/s = 1.4 GeV and also 
leads to a a pole at 445 - z221 MeV. 

As we can see, all these different schemes, which uses the lowest order chiral 
Lagrangian and the implementation of exact unitarity, lead invariably to a a pole 
around 450 — i225 MeV, with small fluctuations. 

Our aim here is to follow these lines and use the unitary scheme to generate 
the 7T7T scalar isoscalar amplitude in the unphysical (virtual) region, and later its 
contribution to the NN interaction. 



2 Lowest order contribution in isoscalar exchange 
in the NN interaction 

We first investigate the lowest-order contribution to the irn interaction in the scalar 
isoscalar channel, which is given by the diagrams of Fig. [1] for pp —>■ pp. We leave 
out the box diagram contribution here and just refer the literature. |1J 

Considering the f isospin dependence of the ttNN coupling, one gets from the 
sum of the diagrams in Fig.[l| the combination, 

^7r 7r°-^7r°7r + 2t 7r 7r 0^ 7r+7r - + 2t 7r + 7r -^ 7r 7r + 4t 7r + 7r -_ >7r + 7r - . (1) 

Taking into account the unitary normalization of the 7T7T states of Ref. [16], 

|7r7r, / = 0) = ^|7T 7T + 7T + 7T^ + 7T-7T+), (2) 

V6 

the combination of t-matrices in Eq. ([!]) corresponds, in terms of the isoscalar tttt 
amplitude, to 

6e = -tr (3) 

We should note here that all the pion lines in Fig. |l| are off shell, which requires the 
use of the off shell t-matrix. From Ref. [16] this off shell amplitude obtained from 
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Figure 1: The lowest-order processes in the 7T7t interaction in the scalar isoscalar 
channel for pp— »pp 



the lowest order meson-meson Lagrangian is given by 

&l = ~ (s* + ^ - 3 £rf) , (4) 

where / is the pion decay constant (/ = 93 MeV), s is the nir Mandelstam variable, 
and the Pi are the momenta of the pion lines. We define the on-shell value of the 
amplitude as that of Eq. (|) when pf = m%. This allows us to write in a convenient 

way 

M=0,L) _ f (I=0,L,OS) , J_ V-/ 2 _ 2n /k\ 

''TTTT^TTTT ~~ 7T7T — >7T7T 3/ 2 . 

with the on shell value of the amplitude 

e^ os) = -i (* - ?fj . (6) 

In what follows we prove that the off shell contribution of the amplitude, i.e., 
the terms proportional to pf — m%, cancels exactly with the diagrams of Fig. |2|, 
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(a) (b) (c) (d) 

Figure 2: The processes with three meson vertex at a baryon line for pp — » pp. 



which appear at the same order of the chiral counting. The diagrams of Fig. [| 
contain the vertices with one baryon line and three mesons. The baryon-meson chiral 
Lagrangian relates this function to that with one meson attached to the baryon line. 
This interaction term is given by 



D -\- F 

:[ B) = — - — {pY^ufp + nYlsufn + nYlsufp + VY^ufi 



where p and n represent the proton or neutron fields. In Eq. (0) is the SU (2) 
matrix given by 



u„ = -y<^$ + ^(d M <M> 2 - 2$^$$ + $ 2 d M $), (8) 

and $ is the £77(2) matrix for the pions given by 



W 71 



*=M£ ■ 0) 
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Equation(^), with the first term of in Eq. (^) gives rise to the ttNN coupling, 



which in the nonrelativistic reduction, 7^75 — > u k 5nki reads 



— tt^NN — t>i — — °" - Pi 

C n + = C n - = \Pl\ C n o pp = 1; C n o nn = — 1 (10) 
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Figure 3: Diagrams involving the three pion vertex. The figure displays the coefli- 

2/ 



cient which multiplies ^j-op for each three pion vertex. 



for incoming pion with momentum p. 

The second term in the expansion of gives rise to the three pion vertices and 
involves derivatives in the three pion fields. However, given the fact that the irNN 
vertex in the loops with meson momentum p of Fig. ^| contain the p-wave coupling 
Bp, only the term which involves the derivative in the pion with momentum p, which 
will give rise to another ap vertex, will contribute in the loop integration. In Fig. [3] 
we show the diagrams corresponding to the loop in Fig. 0(a) for the different isospin 
combinations of Fig. [TJ, together with the contribution of the three meson vertex to 
each diagram. 

If one sums the contributions of all the diagrams in Fig. including the isospin 
weight of the nNN vertices one obtains an equivalent three pion vertex, 

1 D + F^ ^ . 6 D + F^ _ , , 

~ lt -'^W~ a ' p (2 + 2 + 21 = w-— 

Now let q be the momentum exchanged between the protons and let us evaluate 
the diagrams of Fig. [L] with the off shell part of the meson-meson vertex and the 
diagrams of Fig. |3|. We present them in Fig. ^ with the appropriate momentum 
assignment. 

The loops to the right in both diagrams of the figure are the same. Hence, we 
concentrate only on the loop of the left. Let us pick up the term ((p + q) 2 — m 2 )/3f 2 
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Figure 4: Diagrams where off-shell cancellations appear. Diagram (A) with the off 
shell part of the 7r 7r amplitude from the line p + q cancels diagram (B) . 



of the meson-meson vertex, which provides the off shell contribution of the left top 
meson in diagram (A). We then get 



V>A) - 'fir) l^-v + * 8 -'^ + tf-^w=k+* 

(12) 



(p + q) 2 -ml + ie E(p[ + p) p\ + p° - E(p~[ + p) + ie" 

with M and E the mass and on shell energy of the nucleon. Similarly, we obtain 
the contribution from the loop on the left in diagram (B), 

rr(B) .( D + F \ 2 /■ d 4 p _ 1 1 



2f ) J (2vr) 4 * * 3f 2 p 2 - m n 2 + ie 
M 1 



E(p 1 +p) p\ + p° - E(p~[ +p) +ie' 



(13) 



In Eq. di~2[) one can see a cancellation of the off-shell part of the meson-meson vertex 
with the corresponding meson propagator. The rest of the integrand only contains q 
in the term a ■ q a ■ p. This term vanishes exactly in the limit p~[ — > 0, or in any case 
if one makes the heavy baryon approximation pi — E{p[ +p) = 0. (The corrections 
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are of order (pl/M) 2 , which we can neglect.) We shall evaluate the NN potential 
for pi = for simplicity. With the a ■ q a ■ p term of Eq. (12) vanishing, we observe 
that and are equal but with opposite sign, and hence there is an exact 
cancellation of these two terms. 

If now we take the off shell part of the meson vertex corresponding to the other 
mesons we would observe an exact cancellation with the diagrams (b)-(d) of Fig. 
|||. This is interesting for practical purposes since it means at the end that we 
must evaluate only the diagrams of Fig. [I] and using only the on-shell value of the 
meson-meson vertex. 

The need to include this subset of chiral diagrams to find cancellations in the 
NN isoscalar interaction was already stressed in Ref . [9] . There, different arguments 
were used, recalling that in the expansion of the pion field matrix U(jf), the third 
order term has ambiguities and the result cannot depend upon them. The on- 
shell amplitude does not depend on these ambiguities, so here we also prove that, 
after summing the terms discussed above, the results do not depend on unknown 
parameters which would affect the off shell meson-meson amplitude. 

The derivation here is particularly useful for our purposes, since in Ref. 16) it 
was also shown that in the construction of the full meson-meson amplitude, only 
the on-shell part of the meson-meson vertex was needed. This allows us to sum 
immediately the set of diagrams that were included in the unitary Bethe-Salpeter 
approach to the scalar meson amplitude and, hence, in the NN interaction we would 
have the set of diagrams shown in Fig. [5| 
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Figure 5: The nucleon-nucleon interaction in the scalar-isoscalar channel, where the 
mi scattering t-matrix is summed up to all orders in the unitary approach. 
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3 Unitary approach for the interacting pions in 
the isoscalar chnnel 



The diagrams of Fig. [5] are easily summed up. First one must substitute the 
on-shell lowest-order meson-meson amplitude of Eq. (|^) with the Bethe-Salpeter 
amplitude, |I6 

t ('=°) = - J ( S ~ ^ (U) 

w ™ f 2 1 + j- 2 {s - ^)G{s) 1 1 

The function G(s) is the loop function with two pion propagators, 



^)='< j4o_ m 2 I z 77. -2 ~ , , ~ , (15) 



d A q 

(2^f q 2 - m^ 2 + ie (P - qf - + ie 

where P is the total momentum of the two-pion system and P 2 = s. In Ref. [16] 
the integral in Eq. (|15|) was regularized with a cut-off in the CM frame of the two 
mesons. Here it is better to work in the CM of the two nucleons, and hence an 
invariant form for G(s) is preferable. This can be done by using the results for G(s) 
obtained with dimensional regularization, which are shown in Appendix A of Ref. 
[20] to be equivalent to those of a cut-off. We then have 



G(s) 



(4tt) 



-1 + In — — + oTn - 

ji z o — 1 



(16) 



with fi the regularization mass, which was found to be \i = 1.2g max = 1.1 GeV for 
the value of the cutoff, g max , needed for a good fit to the data with the B-S equation 
and the lowest order xPT as the kernel of the equation. In Eq. (0) the magnitude 
of a is 



(17) 

In the range < s < 4m 7r 2 , one has an analytical extrapolation which leads to 

cr+1 /4m 7r 2 /4m 7r 2 

oTn >{ix — 2a)\ 1; a = arctan \ 1, (18) 

a — 1 V s V s 




and for s > Am n 2 the log of Eq. (|16|) develops a negative imaginary part. For s < 0, 
a is always larger than 1, and the log term behaves smoothly. 
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Figure 6: The nucleon-nucleon interaction with the two-pion nucleon vertex in the 
isoscalar channel. 



We would also like to note that at the same order as calculated in this section 
there are other pieces, like in Fig. |6|. The BBMM contact term is proportional to 
7 / V a , which implies vector and iso-vector exchange in the t channel, and hence we 
do not have to consider it for the isoscalar exchange channel we are concerned with 
here. 

In the iViV CM frame we have the momentum q exchanged between the nucleons 
with q = (0, q), such that s = —q 2 . The iViV potential corresponding to the isoscalar 
27r exchange of Fig. [5] then becomes 

^ = ^4 i- G( 'yffW) - (i9> 

Here V(q) is the vertex corresponding to the triangle loop to the left in Fig. |J which 
in the limit p% — > can be written as 



d 4 p (D + F x2 



" x> ' J (27r) 4 \ 2f j (^ +( ^" ^p2_ m ^2j t _^ e _|_ q y _ m _2 _|_ i t 

Ml 

E(fj M + p°-E(p)+ie' ^ 
With the sums over the spins of the intermediate nucleon in the loop we have 
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Figure 7: The two pion exchange triangle vertex through delta excitation. 



a ■ (p + q)a ■ p — (p + q)-p+(qxp)-a, (21) 

but the integral over p in Eq. (|20D leads to a vector proportional to q, which makes 
the term with a of Eq. ([?!]) vanish in the integral. Furthermore, the p° integration 
in Eq. (|20|) can be done analytically, and one finds 



, d 3 p (D + F\ 2 _ 2 _ M 111 1 



(2vr) 3 V 2/ / ^ E 2 w w' w + w' 

' 1 [cuW + £-M], (22) 



£ = E(p); u = y/m^ + p 2 ; J = y/m„ 2 + (p + $ 2 , (23) 



with 



The integral in Eq. (gzj) is logarithmically divergent. 

It is straighforward to include the contribution with an intermediate A instead 
of a nucleon in the vertex function, as depicted in Fig. ^. It can be evaluated in 
a straighforward way, substituting S{Tj for <JiTj with S and T the spin and isospin 
transition operators. 
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The inclusion of the A in the fermionic loops is required in order to satisfy the 
basic rules in the large iV c limit, [21] and we include it here. We obtain 




uj + uj' E A + u- M E A + u' 



E a 2ujuj' 
[uj + uj' + E A - M], (24) 



with M A the A mass and E A = (M A 2 +p 2 ) 1 ^ 2 . We take the empirical value for the 
ratio of the couplings f* NA to f n NN as 2.12 . 

The integral of Eq. (|24]) is also logarithmically divergent. In order to regularize 
the integrals of Eqs. fl2"2"|) and (0), one needs a cutoff or a form factor. Within 
a quark model one finds natural form factors due to the finite size of the nucleon, 
but one must still sum over excited states of the quarks in the intermediate states, 



which leads to a divergence in spite of the form factors. |22| The regularization of 
the loop in this case can be accomplished by a cutoff in the space of intermediate 
states, together with a form factor. The truncation at A, together with the use 
of monopole form factors with A of order 1 GeV has been found to be a sensible 
regularization procedure, as proved by the success of such an approach in the cloudy 
bag model. [23| We shall use here static form factors to keep the p° integral of Eq. 



( |20| ) simple, and hence we include in the integrand of Eqs. ( p2|) and ( |24"D the product 
of form factors 

W«>a4aWF' (25) 

with A » 1 - 1.2 GeV. 

Our final form for the AW potential in momentum space is thus given by 

- ^S0T^) > (26) 

with 

V(q) = V N (q) + V A (q). (27) 
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Here V/v(g) and Va(<?) are given by Eqs. (^) and (|24| ) incorporating in the integrand 
the product F(p)F(p + q) of Eq. (|25|). The potential in coordinate space is given by 

V${r) = J-0^^t NN (Q) 
1 1 f°° 

= / qdqsm(qr)t NN (q). (28) 

We would like to compare vjf^(r) with the empirical a exchange given by 



d 3( l Jqr 9aNN 



NNK ' J ~ / _7?2 _ m 2 



(27r) 3 — q 2 — m 

„2 p -m CT r 

(29) 



with 

n 2 
Act AW 



47T 



5.69; QaNN = 8.46, (30) 



47T 

from Ref. [1]. This corresponds to the empirical a exchange when the two pion box 
diagrams with the intermediate A are considered in addition. 

Actually, in Ref. [1] a monopole form factor per vertex with A a = 1.7 GeV is 
also included so that 



V ia) (a) --a 2 ( A °- m2 A 2 1 ( 3i ) 



which leads to the potential in coordinate space 



V {a) (r) - — fl 2 \-P' h " r 4- ^2_^lp-A CT r _ \-m a r \ foo] 

WJ- 47T 9aNN \r 2A(j r J ' ld j 

with the same asymptotic behaviour as Eq. (|29"D at large distances and the finite 
limit — ( A(j - m a ) 2 /2A a for r -> 0. 
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4 Qualitative discussion of the potential 



Let us go back to Eq. flT4|) in the physical region. We know from Ref. [16] that the 
isoscalar amplitude develops a pole around 450 — z225 MeV. Let us neglect for this 
study ImG(s) and rewrite Eq. ( |i~4"|) as 

.(i=o) _ -G(s)-\s-^) 

so that 



2 i 2 



2 /^(s)- 1 ~ G{s)~' = -('-^f- - m CT 2 ). (34) 

Assuming now that G(s)^ 1 is a smooth function of s and V"(g) is a smooth function 
of g, we would have 



t NN (q) ~ V(0) 



2 6 (m g 2 -^)(g 2 + ^) 

/ 2 (f 2 + m CT 2 

2 6(m CT 2 - ^f) (q 2 + m a 2 ) - (m 2 
/ 2 g 2 + m CT s 



1/(0) " ~ ; W ^ ~ l , "~ - ■ (35) 



We see that the first term in the numerator of the second fraction gives rise to 
a repulsive 5 function in coordinate space, while the second term gives rise to an 
attractive potential due to a standard a exchange with the equivalent coupling 



Tl 2 \ 

g aNN ~ t/(0)V6 Vm ^ ~ ' 2 J (36) 



2 

7 

With the value V"(0) ~ 0.10 x 10~ 2 MeV -1 that we obtain, as we shall see, the 
uNN coupling would be of the order of g a NN — 5, which has the right order of 
magnitude, compared with the empirical coupling used in Ref. [1] to provide the 
needed intermediate range attraction, g a NN — 8.5. 

This qualitative study is certainly improved by the accurate results which we 
give in the next section, but it serves to illustrate the features which we can expect 
from such a potential. The interesting thing is that, apart from the intermediate 
attraction, one gets a scalar repulsion, which is also demanded by the NN scattering 
data. 
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Figure 8: The vertex function V(q) as a function of momentum q. 



In Fig. |8| we give the results for the vertex function, which display a moderate but 
steady decrease of the vertex with q. This ensures convergence of the integral of Eq. 
for V(r). At this point we note already some discrepency with Ref. [9], where 
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□ 



□ 



Figure 9: The function ijvjv- 



( |26|) is plotted in Fig. |9|. It exhibits an increase from to about 500 MeV and then 
drops as q increases further. For q around 2 GeV, it has a value about 250 times 
smaller than at the maximum. 

We should caution that the unitary theory used here leads to good results in 
the physical region up to about y/s = 1.2 GeV. It is then clear that one should not 
extrapolate the results of the model in the unphysical region to values of q much 
larger than 1 GeV. In this sense, the form factors in the nNN vertices with lambda 
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Figure 10: The NN potential in coordinate space 



of order 1 GeV guarantee that one does not enter into this unknown region. In any 
case, the uncertainties for these large values of q would revert into uncertainties at 
short distances in V(r), which we certainly must admit. Technically, the function 
£jV7v(<?) °f Eq. (|26| ) develops a pole around 2060 MeV, far away from the region 
of validity of the model used, and to which we do not give any physical meaning. 
This is anyway a warning that we should not attempt to investigate the behaviour 
of the potential at very short distances with the present model. The numerical 
integrations are extended up to 2 GeV in the actual calculations. The predictions 
for V(r) beyond r = 0.5 fm should be rather safe. 

In Fig. [lO] we show the results for V(r). We find interesting behaviour with a 
moderate attraction beyond r = 0.9 fm and a repulsion at shorter distances, which 
become of order 1 GeV at very short distances, with the caveat discussed above. 
Around 0.5 fm it has a value of around 200 MeV. The trend of the results agrees with 
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the qualitative results anticipated in §4. The first thing to note is that these results 
are quite different qualitatively from those of the conventional sigma exchange, where 
there is only attraction, according to Eq. (P2"|). They also differ appreciably from 
those of Ref. [9] where a very weak repulsion is obtained in the region where we 
find an attraction here (the attraction of the order of eV found numerically at larger 
distances in Ref. 9 is not significative). This finding is relevant since it shows that 
even for large distances, which is the focus of Ref. [9], one needs to go beyond 
perturbation theory in the tttt interaction in order to determine the behaviour of 
the potential. A study of the behaviour of V(r) around 2.5 fm reveals Yukawa 
behaviour with vn a around 450 MeV and a coupling g a NN of about 8, in rough 
agreement with the results of the qualitative discussion of §4, demonstrating that 
the sigma pole appearing in the physical region shows up indeed in the potential 
at large distances. Yet, the presence of the sigma pole is not the only element 
responsible for the behaviour of the potential seen in Fig. [II]. Indeed, if we replace 
the denominator of Eq. (Efj) by unity, thus removing the sigma pole, we still find 
structure of the potential similar to that in Fig. |l(| with the attraction appearing 
at shorter distances beyond r = 0.85 fm. The strength of the attraction is larger 
and also falls off faster with increasing r, indicating that the range is now given by 
other g-dependent functions. In this case, the vertex function and its q dependence 
are responsible for that behaviour and for the differences with the results of Ref. 
[9]. Both the vertex function and the unitary approach to the irn interaction are 
thus very important for the isoscalar potential and for determining its strength in 
all ranges of distances. Note that the phenomenological form factor that we assume 
modulates the structure of the vertex. It not only affects the vertex at large q, 
in which case only the short range part of the interaction would be modified, but 
it also modifies the strength of the vertex at finite q, and this also has influence 
in the medium and long range parts of the interaction. For similar reasons, the 
unitarization changes the q dependence, but more importantly it introduces a pole 
in the s channel, which is reflected in the t channel by a distinct tail in r that goes 
roughly as e~ m " r jr. 

The evaluation of the vertex function here has been done using elements of phe- 
nomenology, beyond the chiral approach otherwise used. In this sense the monopole 
form factors were used and the space of intermediate baryon states was truncated 
in the delta. Other more realistic options could be used, maybe even using form 
factors obtained within the chiral approach in some self-consistent way. For the time 
being we should admit uncertainties from this source. These uncertainties should not 
be minimized since the results depend appreciably on the choice of the range pa- 
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rameters. For instance a change of lambda from 1.1 GeV to 1.2 GeV changes the 
attraction from 10 to 15 MeV. A different vertex evaluation, in which an explicit 
contribution of states beyond the delta in the intermediate states is allowed, as in 
Ref. [22] would also lead to different q behaviour of the vertex function, with imme- 
diate results in the potential in coordinate space. The regularization procedure can 
be done in many ways. In chiral perturbation, it requires a regularization scheme by 
means of a regularization scale or a cut-off, but at the same time one has to intro- 
duce counter-terms of higher order. In the case of the meson-meson interactions, it 
is proved in Ref. [17] and [20] that a cut-off, or equivalently dimensional regulariza- 
tion at a certain scale, effectively generates the counter-terms in the meson-meson 
s-wave channel. For the vertex functions, which we generate here, such a procedure 
is not yet available. In Ref. [8] and [9], infinities in the vertex function are removed 
with the claim that this removal affects only the short range part of the interaction. 
In our case we have resorted to phenomenology using an empirical ttNN form factor 
which makes the loop function convergent. The fact that with this procedure we 
find an assymptotic behavior for the potential similar to the empirical a potential 
gives us confidence that this "empirical" regularization in the baryon sector is a 
reasonable procedure. 

In any case, given all the limitations which have been discussed above, we can 
not make strong claims about the strength of this a potential, but we certainly 
can stress the qualitative features, which are rather different from the conventional 
results for sigma exchange and from results obtained in the perturbative approach 
to the problem. 

6 Conclusions 

We have evaluated the contribution to the NN potential from the exchange of an 
interacting pion pair in the scalar-isoscalar channel, the sigma channel. We have 
used a unitary model based upon the lowest order chiral Lagrangian and the Bethe- 
Salpeter equation. This model is supported by a more elaborate approach using 
the inverse amplitude method and the lowest and second order chiral Lagrangians, 
but technically simpler. Both approaches generate a sigma pole around 450 MeV, 
although with a large width. The tttt amplitude extrapolated to negative values of s, 
retains memory of the sigma pole and leads to a tail of the NN potential reflecting 
a Yukawa behaviour with the range of the sigma mass. We obtain some moderate 
attraction beyond 0.9 fm and a stronger repulsion before that. This behaviour is 
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quite different from that of the conventional sigma exchange, which always leads to 
an attractive force, even in the presence of form factors for the aNN vertices. 

We have also found that the procedure chosen to regularize the aNN vertex 
function also has important consequences for the potential. We have chosen to 
introduce phenomenological form factors for the nNN vertices and to restrict the 
space of baryon intermediate states in the loop to the nucleon and the delta, a 
procedure which has been often used in related problems in connection with chiral 
quark models. The results are sensitive to the form factors and the manner in 
which the vertex is constructed. This is one of the problems that requires further 
thought. For the moment, we should accept one degree of freedom due to the 
divergence of the vertex loop, which must be cured with some regularization mass, 
a counterterm, the use of a form factor, etc. We can hope to reduce the freedom 
by using some empirical information, as in the present case . Hence we believe our 
results are qualitatively acceptable in the range of distance beyond 0.5 fm, and more 
quantitatively accurate results. The findings of this paper are important because 
they reveal behaviour that is quite different from that which one gets from the 
conventional sigma exchange picture. On the other hand, it is also known that a 
nonperturbative approach to the im interaction is needed, and such an approach 
should be important for obtaining proper results in the whole range of distances in 
the NN interaction. These findings should thus have important consequences in the 
chiral approach to the NN interaction, which is now the subject of much theoretical 
investigation. 
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